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Abstract 

It is known that evaluating a certain approximation to the Jones polynomial for the plat closure of 
a braid is a BQP-complete problem. That is, this problem exactly captures the power of the quantum 
circuit model [T3ll3][T|. The one clean qubit model is a model of quantum computation in which all but one 
qubit starts in the maximally mixed state. One clean qubit computers are believed to be strictly weaker 
than standard quantum computers, but still capable of solving some classically intractable problems [21) . 
Here we show that evaluating a certain approximation to the Jones polynomial at a fifth root of unity 
for the trace closure of a braid is a complete problem for the one clean qubit complexity class. That is, 
a one clean qubit computer can approximate these Jones polynomials in time polynomial in both the 
number of strands and number of crossings, and the problem of simulating a one clean qubit computer 
is reducible to approximating the Jones polynomial of the trace closure of a braid. 

1 One Clean Qubit 

The one clean qubit model of quantum computation originated as an idealized model of quantum computation 
on highly mixed initial states, such as appear in NMR implementations [2 1| H] . In this model, one is given an 
initial quantum state consisting of a single qubit in the pure state |0), and n qubits in the maximally mixed 
state. This is described by the density matrix 

P=|0> (0|®|^. 

One can apply any polynomial-size quantum circuit to p, and then measure the first qubit in the com- 
putational basis. Thus, if the quantum circuit implements the unitary transformation U , the probability of 
measuring |0) will be 

po = Tr[(|0) (0| ® I)UpU^] = 2-"Tr[(|0) (0| ® I)U{\Q) (0| ® I)U\ (1) 

Computational complexity classes are typically described using decision problems, that is, problems 
which admit yes/no answers. This is mathematically convenient, and the implications for the complexity of 
non-decision problems are usually straightforward to obtain (c/. [H])- The one clean qubit complexity class 
consists of the decision problems which can be solved in polynomial time by a one clean qubit machine with 
correctness probability of at least 2/3. The experiment described in equation [T] can be repeated polynomially 
many times. Thus, if pi > 1/2 + e for instances to which the answer is yes, and pi < 1/2 — e otherwise, then 
by repeating the experiment poly(l/e) times and taking the majority vote one can achieve 2/3 probability of 
correctness. Thus, as long as e is at least an inverse polynomial in the problem size, the problem is contained 
in the one clean qubit complexity class. Following [21], we will refer to this complexity class as DQCl. 

A number of equivalent definitions of the one clean qubit complexity class can be made. For example, 
changing the pure part of the initial state and the basis in which the final measurement is performed does 
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Figure 1: This circuit implements the Hadamard test. A horizontal line represents a qubit. A horizontal line with 
a slash through it represents a register of multiple qubits. The probability po of measuring |0) is as shown above. 
Thus, one can obtain the real part of {tp\ U lip) to precision e by making 0(l/e^) measurements and counting what 
fraction of the measurement outcomes are |0). Similarly, if the control bit is instead initiaUzed to -^(|0) — i one 
can estimate the imaginary part of {^f}\U {tp). 



not change the resulting complexity class. Less trivially, allowing logarithmically many clean qubits results 
in the same class, as discussed below. It is essential that on a given copy of p, measurements are performed 
only at the end of the computation. Otherwise, one could obtain a pure state by measuring p thus making 
all the qubits "clean" and re-obtaining BQP. Remarkably, it is not necessary to have even one fully polarized 
qubit to obtain the class DQCl. As shown in [21], a single partially polarized qubit suffices. 

In the original definition[21 of DQCl it is assumed that a classical computer generates the quantum 
circuits to be applied to the initial state p. By this definition DQCl automatically contains the complexity 
class P. However, it is also interesting to consider a slightly weaker one clean qubit model, in which the 
classical computer controlling the quantum circuits has only the power of NCI. The resulting complexity 
class appears to have the interesting property that it is incomparable to P. That is, it is not contained in 
P nor does P contain it. We suspect that our algorithm and hardness proof for the Jones polynomial carry 
over straightforwardly to this NCl-controUed one clean qubit model. However, we have not pursued this 
point. 

Any 2" x 2" unitary matrix can be decomposed as a linear combination of ?i-fold tensor products of 
Pauli matrices. As discussed in [TV , the problem of estimating a coefficient in the Pauli decomposition of a 
quantum circuit to polynomial accuracy is a DQCl-complete problem. Estimating the normalized trace of 
a quantum circuit is a special case of this, and it is also DQCl-complete. This point is discussed in ^7\. To 
make our presentation self-contained, we will sketch here a proof that trace estimation is DQCl-complete. 
Technically, we should consider the decision problem of determining whether the trace is greater than a given 
threshold. However, the trace estimation problem is easily reduced to its decision version by the method of 
binary search, so we will henceforth ignore this point. 

First we'll show that trace estimation is contained in DQCl. Suppose we are given a quantum circuit on 
n qubits which consists of polynomially many gates from some finite universal gate set. Given a state \ip) 
of n qubits, there is a standard technique for estimating {ijj\U \ ip), called the Hadamard test [5], as shown in 
figure [T] Now suppose that we use the circuit from figure [l] but choose \ip) uniformly at random from the 
2" computational basis states. Then the probability of getting outcome |0) for a given measurement will be 

1 ^ l + Re{{x\U\x)) _l Re(Tr [/) 
xefo,!}" 

Choosing {ip) uniformly at random from the 2" computational basis states is exactly the same as inputting 
the density matrix //2" to this register. Thus, the only clean qubit is the control qubit. Trace estimation 
is therefore achieved in the one clean qubit model by converting the given circuit for U into a circuit for 
controlled-?/ and adding Hadamard gates on the control bit. One can convert a circuit for U into a circuit 
for controUed-C/ by replacing each gate G with a circuit for controUed-G. The overhead incurred is thus 
bounded by a constant factor [23]. 

Next we'll show that trace estimation is hard for DQCl. Suppose we are given a classical description of 
a quantum circuit implementing some unitary transformation U on n qubits. As shown in equation [l] the 
probability of obtaining outcome |0) from the one clean qubit computation of this circuit is proportional to 
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Figure 2: Here CNOT gates are used to simulate 3 clean ancilla qubits. 



the trace of the non-unitary operator (|0) (0| (E) I)U{\0) (0| ^ I)W , which acts on n qubits. Estimating this 
can be achieved by estimating the trace of 
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which is a unitary operator on n + 2 qubits. This suffices because 

Tr[(|0) (0| ® I)Um (0| ® /)[/t] = -Tt[U']. (2) 
To see this, we can think in terms of the computational basis: 

Tr[[/']= W\x). 

a:G{0,l}" 

If the first qubit of \x) is |I), then the rightmost CNOT in U' will flip the lowermost qubit. The resulting 
state will be orthogonal to \x) and the corresponding matrix element will not contribute to the trace. Thus 
this CNOT gate simulates the initial projector |0) (0| ® / in equation [51 Similarly, the other CNOT in U' 
simulates the other projector in equation [21 

The preceding analysis shows that, given a description of a quantum circuit implementing a unitary 
transformation U on n-qubits, the problem of approximating ^Tr U to within =fc poiy(„) precision is DQCl- 
complete. 

Some unitaries may only be efficiently implementable using ancilla bits. That is, to implement U on 
n-qubits using a quantum circuit, it may be most efficient to construct a circuit on n + m qubits which acts 
as [/ (8) /, provided that the m ancilla qubits are all initialized to |0). These ancilla qubits are used as work 
bits in intermediate steps of the computation. To estimate the trace of [/, one can construct a circuit Ua on 
n + 2m qubits by adding CNOT gates controlled by the m ancilla qubits and acting on m extra qubits, as 
shown in figure [21 This simulates the presence of m clean ancilla qubits, because if any of the ancilla qubits 
is in the |1) state then the CNOT gate will flip the corresponding extra qubit, resulting in an orthogonal 
state which will not contribute to the trace. 

With one clean qubit, one can estimate the trace of Ua to a precision of poiy(„ ^) ■ I^y construction, 

Tr[[/Q] = 2™Tr[[/]. Thus, if m is logarithmic in n, then one can obtain Tr[i7] to precision pofy(„-) , just as can 
be obtained for circuits not requiring ancilla qubits. This line of reasoning also shows that the fc-clean qubit 
model gives rise to the same complexity class as the one clean qubit model, for any constant fc, and even for 
k growing logarithmically with n. 

It seems unlikely that the trace of these exponentially large unitary matrices can be estimated to this 
precision on a classical computer in polynomial time. Thus it seems unlikely that DQCl is contained in 
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Figure 3: Shown from left to right are the unknot, another representation of the unknot, an oriented trefoil knot, 
and the Hopf link. Broken lines indicate undercrossings. 



Figure 4: Two knots are the same if and only if one can be deformed into the other by some sequence of the three 
Reidemeister moves shown above. 



P. (For more detailed analysis of this point see [lO].) However, it also seems unlikely that DQCl contains 
all of BQP. In other words, one clean qubit computers seem to provide exponential speedup over classical 
computation for some problems despite being strictly weaker than standard quantum computers. 

2 Jones Polynomials 

A knot is defined to be an embedding of the circle in M"^ considered up to continuous transformation (isotopy) . 
More generally, a link is an embedding of one or more circles in M.^ up to isotopy. In an oriented knot or 
link, one of the two possible traversal directions is chosen for each circle. Some examples of knots and links 
are shown in figure [3l One of the fundamental tasks in knot theory is, given two representations of knots, 
which may appear superficially different, determine whether these both represent the same knot. In other 
words, determine whether one knot can be deformed into the other without ever cutting the strand. 

Reidemeister showed in 1927 that two knots are the same if and only if one can be deformed into the 
other by some sequence constructed from three elementary moves, known as the Reidemeister moves, shown 
in figured! This reduces the problem of distinguishing knots to a combinatorial problem, although one for 
which no efficient solution is known. In some cases, the sequence of Reidemeister moves needed to show 
equivalence of two knots involves intermediate steps that increase the number of crossings. Thus, it is 
very difficult to show upper bounds on the number of moves necessary. The most thoroughly studied knot 
equivalence problem is the problem of deciding whether a given knot is equivalent to the unknot. Even 
showing the decidability of this problem is highly nontrivial. This was achieved by Haken in 1961 [Mj. In 
1998 it was shown by Hass, Lagarias, and Pippenger that the problem of recognizing the unknot is contained 



A knot invariant is any function on knots which is invariant under the Reidemeister moves. Thus, a 
knot invariant always takes the same value for different representations of the same knot, such as the two 
representations of the unknot shown in figureO In general, there can be distinct knots which a knot invariant 
fails to distinguish. 

One of the best known knot invariants is the Jones polynomial, discovered in 1985 by Vaughan Jones [18] . 
To any oriented knot or link, it associates a Laurent polynomial in the variable t^^"^. The Jones polynomial 
has a degree in t which grows at most linearly with the number of crossings in the link. The coefficients 
are all integers, but they may be exponentially large. Exact evaluation of Jones polynomials at all but a 
few special values of t is #P-hard[T6]. The Jones polynomial can be defined recursively by a simple "skein" 
relation. However, for our purposes it will be more convenient to use a definition in terms of a representation 
of the braid group, as discussed below. 

To describe in more detail the computation of Jones polynomials we must specify how the knot will 




in NP[T5]. 
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Figure 5: Shown from left to right are a braid, its plat closure, and its trace closure. 



be represented on the computer. Although an embedding of a circle in MP is a continuous object, all the 
topologically relevant information about a knot can be described in the discrete language of the braid group. 
Links can be constructed from braids by joining the free ends. Two ways of doing this are taking the 
plat closure and the trace closure, as shown in figure [51 Alexander's theorem states that any link can be 
constructed as the trace closure of some braid. Any link can also be constructed as the plat closure of some 
braid. This can be easily proven as a corollary to Alexander's theorem, as shown in figure [51 





Figure 6: A trace closure of a braid on n strands can be converted to a plat closure of a braid on 2n strands by 
moving the "return" strands into the braid. 
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Figure 7: Shown are the two Markov moves. Here the boxes represent arbitrary braids. If a function on braids is 
invariant under these two moves, then the corresponding function on links induced by the trace closure is a link 
invariant. 



Given that the trace closure provides a correspondence between links and braids, one may attempt to 
find functions on braids which yield link invariants via this correspondence. Markov's theorem shows that 
a function on braids will yield a knot invariant provided it is invariant under the two Markov moves, shown 
in figure [71 Thus the Markov moves provide an analogue for braids of the Reidemeister moves on links. The 
constraints imposed by invariance under the Reidemeister moves are enforced in the braid picture jointly by 
invariance under Markov moves and by the defining relations of the braid group. 
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A linear function / satisfying f{AB) = f{BA) is called a trace. The ordinary trace on matrices is one 
such function. Taking a trace of a representation of the braid group yields a function on braids which is 
invariant under Markov move I. If the trace and representation are such that the resulting function is also 
invariant under Markov move II, then a link invariant will result. The Jones polynomial can be obtained in 
this way. 

In [5], Aharonov, et al. show that an additive approximation to the Jones polynomial of the plat or trace 
closure of a braid aX t = e*^'^/'^ can be computed on a quantum computer in time which scales polynomially 
in the number of strands and crossings in the braid and in k. In [U [3T], it is shown that for plat closures, 
this problem is BQP-complete. The complexity of approximating the Jones polynomial for trace closures 
was left open, other than showing that it is contained in BQP. 

The results of (31 [TJ [JT] reformulate and generalize the previous results of Freedman et al. [T51 [T^ . 
which show that certain approximations of Jones polynomials are BQP-complete. The work of Freedman 
et al. in turn builds upon Witten's discovery of a connection between Jones polynomials and topological 
quantum field theory [30] . Recently, Aharonov et al. have generalized further, obtaining an efficient quantum 
algorithm for approximating the Tutte polynomial for any planar graph, at any point in the complex plane, 
and also showing BQP-hardness at some points [5]. As special cases, the Tutte polynomial includes the 
Jones polynomial, other knot invariants such as the HOMFLY polynomial, and partition functions for some 
physical models such as the Potts model. 

The algorithm of Aharonov et al. works by obtaining the Jones polynomial as a trace of the path model 
representation of the braid group. The path model representation is unitary for t = e*^'^/'^ and, as shown 
in [3], can be efficiently implemented by quantum circuits. For computing the trace closure of a braid the 
necessary trace is similar to the ordinary matrix trace except that only a subset of the diagonal elements of 
the unitary implemented by the quantum circuit are summed, and there is an additional weighting factor. For 
the plat closure of a braid the computation instead reduces to evaluating a particular matrix element of the 
quantum circuit. Aharonov et al. also use the path model representation in their proof of BQP-completeness. 

Given a braid &, we know that the problem of approximating the Jones polynomial of its plat closure 
is BQP-hard. By Alexander's theorem, one can obtain a braid h' whose trace closure is the same link as 
the plat closure of h. The Jones polynomial depends only on the link, and not on the braid it was derived 
from. Thus, one may ask why this doesn't immediately imply that estimating the Jones polynomial of the 
trace closure is a BQP-hard problem. The answer lies in the degree of approximation. As discussed in 
section ini the BQP-complete problem for plat closures is to approximate the Jones polynomial to a certain 
precision which depends exponentially on the number of strands in the braid. The number of strands in b' 
can be larger than the number of strands in 6, hence the degree of approximation obtained after applying 
Alexander's theorem may be too poor to solve the original BQP-hard problem. 

The fact that computing the Jones polynomial of the trace closure of a braid can be reduced to estimating 
a generalized trace of a unitary operator and the fact that trace estimation is DQCl-complete suggest a 
connection between Jones polynomials and the one clean qubit model. Here we find such a connection by 
showing that evaluating a certain approximation to the Jones polynomial of the trace closure of a braid 
at a fifth root of unity is DQCl-complete. The main technical difficulty is obtaining the Jones polynomial 
as a trace over the entire Hilbert space rather than as a summation of some subset of the diagonal matrix 
elements. To do this we will not use the path model representation of the braid group, but rather the 
Fibonacci representation, as described in the next section. 

3 Fibonacci Representation 

The Fibonacci representation p^p^ of the braid group Bn is described in [19] in the context of Temperley-Lieb 
recoupling theory. Temperley-Lieb recoupling theory describes two species of idealized "particles" denoted 
by p and We will not delve into the conceptual and mathematical underpinnings of Temperley-Lieb 
recoupling theory. For present purposes, it will be sufficient to regard it as a formal procedure for obtaining 
a particular unitary representation of the braid group whose trace yields the Jones polynomial a,tt = e*^'^/^. 
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Figure 8: For an n-strand braid we can write a length n + 1 string of p and * symbols across the base. The string 
may have no two * symbols in a row, but can be otherwise arbitrary. 
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Figure 9: ai denotes the elementary crossing of strands i and i + The braid group on n strands B„ is generated 
by o"! . . . (T„_i, which satisfy the relations UifJj = (Tja for \i — j\ > 1 and Gij^iaiOi^-i — UiUi^iUi for all i. The group 
operation corresponds to concatenation of braids. 



Throughout most of this paper it will be more convenient to express the Jones polynomial in terms of 
A = e-^S'^/s, with t defined by t = A-^. 

It is worth noting that the Fibonacci representation is a special case of the path model representation 
used in [3]. The path model representation applies when t = e*^^/*-' for any integer k, whereas the Fibonacci 
representation is for fc = 5. The relationship between these two representations is briefly discussed in 
appendix [C] However, for the sake of making our discussion self contained, we will derive all of our results 
directly within the Fibonacci representation. 

Given an n-strand braid b € Bn, we can write a length ri + 1 string of p and * symbols across the base 
as shown in figure [S] These strings have the restriction that no two * symbols can be adjacent. The number 
of such strings is /n+3, where /„ is the n}^ Fibonacci number, defined so that /i = 1, /2 = 1, /a = 2, . . . 
Thus the formal linear combinations of such strings form an /„+3-dimensional vector space. For each n, the 

(n) 

Fibonacci representation pp is a homomorphism from _B„ to the group of unitary linear transformations on 
this space. We will describe the Fibonacci representation in terms of its action on the elementary crossings 
which generate the braid group, as shown in figure [9l 

The elementary crossings correspond to linear operations which mix only those strings which differ by 
the symbol beneath the crossing. The linear transformations have a local structure, so that the coefficients 
for the symbol beneath the crossing to be changed or unchanged depend only on that symbol and its two 
neighbors. For example, using the notation of [19], 



P 




+ d 

P * p p p p 



which means that the elementary crossing cr,; corresponds to a linear transformation which takes any string 
whose i'^ through (i + 2)**^ symbols are p * p to the coefficient c times the same string plus the coefficient 
d times the same string with the * at the (i + 1)**^ position replaced by p. (As shown in figure 9, the i^^ 
crossing is over the (i + 1)'^ symbol.) To compute the linear transformation that the representation of a 
given braid applies to a given string of symbols, one can write the symbols across the base of the braid, and 
then apply rules of the form[3]until all the crossings are removed, and all that remains are various coefficients 
for different strings to be written across the base of a set of straight strands. 

For compactness, we will use (p^p) — c{p ^ p) + d{ppp) as a shorthand for equation O In this notation. 
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the complete set of rules is as follows. 



where 



{*PP) 

{P*P) 
(PP*) 
(PPP) 



a 
b 
c 
d 
e 
A 



a{*pp) 

c{p *p) + d{ppp) 
a(jpp*) 

d{p*p) + e{ppp), 



-A" 
A^ 

^8^3/2 ^ ^4^3/2 

T - / 
37r/5 



(4) 



A\ - 



2/(1 + x/5). 



(5) 



Using these rules we can calculate any matrix from the Fibonacci representation of the braid group. 
Notice that this is a reducible representation. These rules do not allow the rightmost symbol or leftmost 
symbol of the string to change. Thus the vector space decomposes into four invariant subspaces, namely the 
subspace spanned by strings which begin and end with p, and the p...*, and * . . .p subspaces. As 

an example, we can use the above rules to compute the action of By, on the * . . .p subspace. 



(3)/ \ 



*p*p 
*PPP 



(3)/ ^ 

P*p (^2) 



In appendix |A] we prove that the Jones polynomial evaluated at t 



d 

e 

„i27r/5 



>l=p>l=p 
*PPP 



(6) 



can be obtained as a weighted 



trace of the Fibonacci representation over the * ... * and * . . .p subspaces. 



4 Computing the Jones Polynomial in DQCl 

As mentioned previously, the Fibonacci representation acts on the vector space of formal linear combinations 
of strings of p and * symbols in which no two * symbols are adjacent. The set of length n strings of this 
type, P„, has fn+2 elements, where /„ is the n"^ Fibonacci number: /i = 1, /2 = 1, /3 = 2, and so on. As 
shown in appendix |E1 one can construct a bijective correspondence between these strings and the integers 
from to ,fn+2 — 1 as follows. If we think of * as 1 and p as 0, then with a string s„s„_i . . . si we associate 
the integer 

n 

zis) =^Sif^+i. (7) 

i=l 

This is known as the Zeckendorf representation. 

Representing integers as bitstrings by the usual method of place value, we thus have a correspondence 
between the elements of P„ and the bitstrings of length b — [log2(/n+2)l • This correspondence will be a 
key element in computing the Jones polynomial with a one clean qubit machine. Using a one clean qubit 
machine, one can compute the trace of a unitary over the entire Hilbert space of 2" bitstrings. Using CNOT 
gates as above, one can also compute with polynomial overhead the trace over a subspace whose dimension 
is a polynomially large fraction of the dimension of the entire Hilbert space. However, it is probably not 
possible in general for a one clean qubit computer to compute the trace over subspaces whose dimension is 
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an exponentially small fraction of the dimension of the total Hilbert space. For this reason, directly mapping 
the strings of p and * symbols to strings of 1 and will not work. In contrast, the correspondence described 
in equation [7] maps P„ to a subspace whose dimension is at least half the dimension of the full 2''-dimensional 
Hilbert space. 

In outline, the DQCl algorithm for computing the Jones polynomial works as follows. Using the results 
described in section [H we will think of the quantum circuit as acting on b maximally mixed qubits plus 
0(1) clean qubits. Thinking in terms of the computational basis, we can say that the first b qubits are in 
a uniform probabilistic mixture of the 2^ classical bitstring states. By equation [71 most of these bitstrings 
correspond to elements of P„. In the Fibonacci representation, an elementary crossing on strands i and z — 1 
corresponds to a linear transformation which can only change the value of the i^^ symbol in the string of p's 
and *'s. The coefficients for changing this symbol or leaving it fixed depend only on the two neighboring 
symbols. Thus, to simulate this linear transformation, we will use a quantum circuit which extracts the 
{i — i*"^, and (z + 1)"^ symbols from their bitstring encoding, writes them into an ancilla register while 
erasing them from the bitstring encoding, performs the unitary transformation prescribed by equation [4] on 
the ancillas, and then transfers this symbol back into the bitstring encoding while erasing it from the ancilla 
register. Constructing one such circuit for each crossing, multiplying them together, and performing DQCl 
trace-estimation yields an approximation to the Jones polynomial. 

Performing the linear transformation demanded by equation [4] on the ancilla register can be done easily 
by invoking gate set universality (c/. Solovay-Kitaev theorem [53]) since it is just a three-qubit unitary 
operation. The harder steps are transferring the symbol values from the bitstring encoding to the ancilla 
register and back. 

It may be difficult to extract an arbitrary symbol from the bitstring encoding. However, it is relatively 
easy to extract the leftmost "most significant" symbol, which determines whether the Fibonacci number /„ 
is present in the sum shown in equation [7] This is because, for a string s of length n, z{s) > fn-i if and 
only if the leftmost symbol is *. Thus, starting with a clean |0) ancilla qubit, one can transfer the value 
of the leftmost symbol into the ancilla as follows. First, check whether z{s) (as represented by a bitstring 
using place value) is > fn-i- If so flip the ancilla qubit. Then, conditioned on the value of the ancilla qubit, 
subtract fn-i from the bitstring. (The subtraction will be done modulo 2^ for reversibility.) 

Any classical circuit can be made reversible with only constant overhead. It thus corresponds to a 
unitary matrix which permutes the computational basis. This is the standard way of implementing classical 
algorithms on a quantum computer [23j. However, the resulting reversible circuit may require clean ancilla 
qubits as work space in order to be implemented efficiently. For a reversible circuit to be implementable 
on a one clean qubit computer, it must be efficiently implementable using at most logarithmically many 
clean ancillas. Fortunately, the basic operations of arithmetic and comparison for integers can all be done 
classically by NCI circuits [29]. NCI is the complexity class for problems solvable by classical circuits of 
logarithmic depth. As shown in [4|, any classical NCI circuit can be converted into a reversible circuit using 
only three clean ancillas. This is a consequence of Barrington's theorem. Thus, the process described above 
for extracting the leftmost symbol can be done efficiently in DQCl. 

More specifically, Krapchenko's algorithm for adding two n-bit numbers has depth [logn] + 0{\/\ogn) 
|29j . A lower bound of depth logn is also known, so this is essentially optimal :29 . Barrington's construction 
[7] yields a sequence of 2^^^ gates on 3 clean ancilla qubits [J] to simulate a circuit of depth d. Thus we obtain 
an addition circuit which has quadratic size (up to a subpolynomial factor). Subtraction can be obtained 
analogously, and one can determine whether a > b can be done by subtracting a from b and looking at 
whether the result is negative. 

Although the construction based on Barrington's theorem has polynomial overhead and is thus sufficient 
for our purposes, it seems worth noting that it is possible to achieve better efficiency. As shown by Draper 
[llj . there exist ancilla-free quantum circuits for performing addition and subtraction, which succeed with 
high probability and have nearly linear size. Specifically, one can add or subtract a hardcoded number 
a into an n-qubit register \x) modulo 2" by performing quantum Fourier transform, followed by O(n^) 
controUed-rotations, followed by an inverse quantum Fourier transform. Furthermore, using approximate 
quantum Fourier transforms [3 [6], [11] describes an approximate version of the circuit, which, for any value 
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Figure 10: Here we make a correspondence between strings of p and * symbols and ordered pairs of integers. The 
string of 9 symbols is split into substrings of length 4 and 5, and each one is used to compute an integer by adding 
the {i + 1)*'^ Fibonacci number if * appears in the place. Note the two strings are read in different directions. 

of parameter m, uses a total of only O^mnlogn) gatef0 to produce an output having an inner product with 
\x + a mod 2") of 1 - 0(2^"). 

Because they operate modulo 2", Draper's quantum circuits for addition and subtraction do not imme- 
diately yield fast ancilla-free quantum circuits for comparison, unlike the classical case. Instead, start with 
an n-bit number x and then introduce a single clean ancilla qubit initialized to |0). Then subtract an n-bit 
hardcoded number a from this register modulo 2"+"'^. If a > a; then the result will wrap around into the 
range [2", 2"+^ — 1], in which case the leading bit will be 1. If a < a; then the result will be in the range 
[0,2" — 1]. After copying the result of this leading qubit and uncomputing the subtraction, the comparison 
is complete. Alternatively, one could use the linear size quantum comparison circuit devised by Takahashi 
and Kunihiro, which uses n uninitialized ancillas but no clean ancillas[28]. 

Unfortunately, most crossings in a given braid will not be acting on the leftmost strand. However, we 
can reduce the problem of extracting a general symbol to the problem of extracting the leftmost symbol. 
Rather than using equation [7] to make a correspondence between a string from P„ and a single integer, 
we can split the string at some chosen point, and use equation [7] on each piece to make a correspondence 
between elements of P„ and ordered pairs of integers, as shown in figure [TUl To extract the i"^ symbol, we 
thus convert encoding [7] to the encoding where the string is split between the i**^ and (i — 1)**^ symbols, so 
that one only needs to extract the leftmost symbol of the second string. Like equation [71 this is also an 
efficient encoding, in which the encoded bitstrings form a large fraction of all possible bitstrings. 

To convert encoding [7] to a split encoding with the split at an arbitrary point, we can move the split 
rightward by one symbol at a time. To introduce a split between the leftmost and second-to-leftmost symbols, 
one must extract the leftmost symbol as described above. To move the split one symbol to the right, one 
must extract the leftmost symbol from the right string, and if it is * then add the corresponding Fibonacci 
number to the left string. This is again a procedure of addition, subtraction, and comparison of integers. 
Note that the computation of Fibonacci numbers in NCI is not necessary, as these can be hardcoded into 
the circuits. Moving the split back to the left works analogously. As crossings of different pairs of strands 
are being simulated, the split is moved to the place that it is needed. At the end it is moved all the way 
leftward and eliminated, leaving a superposition of bitstrings in the original encoding, which have the correct 
coefficients determined by the Fibonacci representation of the given braid. 

Lastly, we must consider the weighting in the trace, as described by equation [TOl Instead of weight Ws, 
we will use Wsl4> so that the possible weights are 1 and 1/0 both of which are < 1. We can impose any 
weight < 1 by doing a controlled rotation on an extra qubit. The CNOT trick for simulating a clean qubit 
which was described in section [T] can be viewed as a special case of this. All strings in which that qubit takes 
the value |1) have weight zero, as imposed by a 7r/2 rotation on the extra qubit. Because none of the weights 
are smaller than 1/0, the weighting will cause only a constant overhead in the number of measurements 
needed to get a given precision. 

5 DQCl-hardness of Jones Polynomials 

We will prove DQCl-hardness of the problem of estimating the Jones polynomial of the trace closure of a 
braid by a reduction from the problem of estimating the trace of a quantum circuit. To do this, we will 
specify an encoding, that is, a map rj : Qn — > Sm from the set Qn of strings of p and * symbols which 

'^A linear-size quantum circuit for exact ancilla-free addition is known, but it does not generalize easily to the case of 
hardcoded summands [9]. 
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start with * and have no two * symbols in a row, to Sm, the set of bitstrings of length m. For a given 
quantum circuit, we will construct a braid whose Fibonacci representation implements the corresponding 
unitary transformation on the encoded bits. The Jones polynomial of the trace closure of this braid, which 
is the trace of this representation, will equal the trace of the encoded quantum circuit. 

Unlike in section [H we will not use a one to one encoding between bit strings and strings of p and * 
symbols. All we require is that a sum over all strings of p and * symbols corresponds to a sum over bitstrings 
in which each bitstring appears an equal number of times. Equivalently, all bitstrings b € Sm must have a 
preimage r]~^(b) of the same size. This insures an unbiased trace in which no bitstrings are overweighted. 
To achieve this we can divide the symbol string into blocks of three symbols and use the encoding 

PPP - 

p*p ^1 ^ ' 

The strings other than ppp and p*p do not correspond to any bit value. Since both the encoded 1 and the 
encoded begin and end with p, they can be preceded and followed by any allowable string. Thus, changing 
an encoded 1 to an encoded zero does not change the number of allowed strings of p and * consistent with 
that encoded bitstring. Thus the condition that be independent of b is satisfied. 

We would also like to know a priori where in the string of p and * symbols a given bit is encoded. This 
way, when we need to simulate a gate acting on a given bit, we would know which strands the corresponding 
braid should act on. If we were to simply divide our string of symbols into blocks of three and write down 
the corresponding bit string (skipping every block which is not in one of the two coding states ppp and p*p) 
then this would not be the case. Thus, to encode n bits, we will instead divide the string of symbols into 
n superblocks, each consisting of clogn blocks of three for some constant c. To decode a superblock, scan 
it from left to right until you reach either a ppp block or a p>i=p block. The first such block encountered 
determines whether the superblock encodes a 1 or a 0, according to equation [8] Now imagine we choose a 
string randomly from Qscniogn- By choosing the constant prefactor c in our superblock size we can ensure 
that in the entire string of Sen log n symbols, the probability of there being any noncoding superblock which 
contains neither a ppp block nor a p>i=p block is polynomially small. If this is the case, then these noncoding 
strings will contribute only a polynomially small additive error to the estimate of the circuit trace, on par 
with the other sources of error. 

The gate set consisting of the CNOT, Hadamard, and tt/8 gates is known to be universal for BQP [23] . 
Thus, it suffices to consider the simulation of 1-qubit and 2-qubit gates. Furthermore, it is sufficient to 
imagine the qubits arranged on a line and to allow 2-qubit gates to act only on neighboring qubits. This 
is because qubits can always be brought into neighboring positions by applying a series of SWAP gates to 
nearest neighbors. By our encoding a unitary gate applied to qubits i and i + 1 will correspond to a unitary 
transformation on symbols iSclogn through (i + 2)3clogn — 1. The essence of our reduction is to take each 
quantum gate and represent it by a corresponding braid on logarithmically many symbols whose Fibonacci 
representation performs that gate on the encoded qubits. 

Let's first consider the problem of simulating a gate on the first pair of qubits, which are encoded in 
the leftmost two superblocks of the symbol string. We'll subsequently consider the more difficult case of 
operating on an arbitrary pair of neighboring encoded qubits. As mentioned in section |3l the Fibonacci 

(n) (n) 

representation pp is reducible. Let pU denote the representation of the braid group _B„ defined by the 
action of on the vector space spanned by strings which begin and end with As shown in appendix 

IbI p)!l\Bn) taken modulo phase is a dense subgroup of SU{fn-i), and p^J^{Bn) modulo phase is a dense 
subgroup of SU{fn). 

In addition to being dense, the ** and *p blocks of the Fibonacci representation can be controlled 
independently. This is a consequence of the decoupling lemma, as discussed in appendix |B] Thus, given 
a string of symbols beginning with *, and any desired pair of unitaries on the corresponding *p and 
vector spaces, a braid can be constructed whose Fibonacci representation approximates these unitaries to 
any desired level of precision. However, the number of crossings necessary may in general be large. The 
space spanned by strings of logarithmically many symbols has only polynomial dimension. Thus, one might 
guess that the braid needed to approximate a given pair of unitaries on the *p and vector spaces for 
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Figure 11: This sequence of unitary steps is used to bring a * symbol where it is needed in the symbol string to 
ensure density of the braid group representation. The presence of the left * ensures density to allow the movement 
of the right * by proposition [1] Similarly, the presence of the right * allows the left * to be moved. 



logarithmically many symbols will have only polynomially many crossings. It turns out that this guess is 
correct, as we state formally below. 

Proposition 1. Given any pair of elements U^,p G SU{fk+i) and U^^, G SU{fk), and any real parameter e, 
one can in polynomial time find a braid b £ Bk with poly(n, log(l/e)) crossings whose Fibonacci representa- 
tion satisfies \\p*p{b) — U^,p\\ < e and \\p^,^,{b) — C/*»|| < e, provided that k — O(logn). By symmetry, the same 
holds when considering pp^, rather than p^,p. 

Note that proposition [T] is a property of the Fibonacci representation, not a generic consequence of 
density, since it is in principle possible for the images of group generators in a dense representation to lie 
exponentially close to some subgroup of the corresponding unitary group. We prove this proposition in 
appendix |D] 

With proposition [1] in hand, it is apparent that any unitary gate on the first two encoded bits can be 
efficiently performed. To similarly simulate gates on arbitrary pairs of neighboring encoded qubits, we will 
need some way to unitarily bring a * symbol to a known location within logarithmic distance of the relevant 
encoded qubits. This way, we ensure that we are acting in the *p or subspaces. 

To move * symbols to known locations we'll use an "inchworm" structure which brings a pair of * symbols 
rightward to where they are needed. Specifically, suppose we have a pair of superblocks which each have a 

* in their exact center. The presence of the left * and the density of p*p allow us to use proposition [T] to 
unitarily move the right * one superblock to the right by adding polynomially many crossings to the braid. 
Then, the presence of the right * and the density of pp^ allow us to similarly move the left * one superblock 
to the right, thus bringing it into the superblock adjacent to the one which contains the right This is 
illustrated in figure [TlJ To move the inchworm to the left we use the inverse operation. 

To simulate a given gate, one first uses the previously described procedure to make the inchworm crawl 
to the superblocks just to the left of the superblocks which encode the qubits on which the gate acts. Then, 
by the density of p^p and proposition [1] the desired gate can be simulated using polynomially many braid 
crossings. 

To get this process started, the leftmost two superblocks must each contain a * at their center. This occurs 
with constant probability. The strings in which this is not the case can be prevented from contributing to 
the trace by a technique analogous to that used in section [T] to simulate logarithmically many clean ancillas. 
Namely, an extra encoded qubit can be conditionally fiipped if the first two superblocks do not both have 

* symbols at their center. This can always be done using proposition [1] since the leftmost symbol in the 
string is always and the p*p and p^^ representations are both dense. 
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swap 



*PPP(5)P*P* >[:pp*(p)*ppp 

Figure 12: This unitary procedure starts with a * in the current superblock and brings it to the center of the target 
superblock. 



It remains to specify the exact unitary operations which move the inchworm. Suppose we have a current 
superblock and a target superblock. The current superblock contains a * in its center, and the target 
superblock is the next superblock to the right or left. We wish to move the * to the center of the target 
superblock. To do this, we can select the smallest segment around the center such that in each of these 
superblocks, the segment is bordered on its left and right by p symbols. This segment can then be swapped, 
as shown in figure [T^ 

For some possible strings this procedure will not be well defined. Specifically there may not be any 
segment which contains the center and which is bordered by p symbols in both superblocks. On such 
strings we define the operation to act as the identity. For random strings, the probability of this decreases 
exponentially with the superblock size. Thus, by choosing c sufhciently large we can make this negligible for 
the entire computation. 

As the inchworm moves rightward, it leaves behind a trail. Due to the swapping, the superblocks are not 
in their original state after the inchworm has passed. However, because the operations are unitary, when the 
inchworm moves back to the left, the modifications to the superblocks get undone. Thus the inchworm can 
shuttle back and forth, moving where it is needed to simulate each gate, always stopping just to the left of 
the superblocks corresponding to the encoded qubits. 

The only remaining detail to consider is that the trace appearing in the Jones polynomial is weighted 
depending on whether the last symbol is p or *, whereas the DQCl-complete trace estimation problem is 
for completely unweighted traces. This problem is easily solved. Just introduce a single extra superblock at 
the end of the string. After bringing the inchworm adjacent to the last superblock, apply a unitary which 
performs a conditional rotation on the qubit encoded by this superblock. The rotation will be by an angle 
so that the inner product of the rotated qubit with its original state where (j) is the golden ratio. This 

will be done only if the last symbol is p. This exactly cancels out the weighting which appears in the formula 
for the Jones polynomial, as described in appendix \K\ 

Thus, for appropriate e, approximating the Jones polynomial of the trace closure of a braid to within ±e 
is DQCl-hard. 



6 Conclusion 

The preceding sections show that the problem of approximating the Jones polynomial of the trace closure of a 
braid with n strands and m crossings to within ±e at t = e'^'^^^ is a DQCl-complete problem for appropriate 
e. The proofs are based on the problem of evaluating the Markov trace of the Fibonacci representation 
of a braid to p^jy^^ precision. By equation 1111 we see that this corresponds to evaluating the Jones 

polynomial with ^ po{^(-„ ,|^) precision, where D = —A^ — A^^ = 2cos(67r/5). Whereas approximating the 
Jones polynomial of the plat closure of a braid was known[l to be BQP-complete, it was previously only 
known that the problem of approximating the Jones polynomial of the trace closure of a braid was in BQP. 
Understanding the complexity of approximating the Jones polynomial of the trace closure of a braid to 

precision ^-^^^jj^-^ was posed as an open problem in [3]. This paper shows that for A = e"*^'^/^, this 
problem is DQCl-complete. Such a completeness result improves our understanding of both the difficulty of 
the Jones polynomial problem and the power one clean qubit computers by finding an equivalence between 
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the two. 

It is generally believed that DQCl is not contained in P and does not contain all of BQP. The DQCl- 
completeness result shows that if this belief is true, it implies that approximating the Jones polynomial of 
the trace closure of a braid is not so easy that it can be done classically in polynomial time, but is not so 
difficult as to be BQP-hard. 

To our knowledge, the problem of approximating the Jones polynomial of the trace closure of a braid is 
one of only four known candidates for classically intractable problems solvable on a one clean qubit computer. 
The others are estimating the Pauli decomposition of the unitary matrix corresponding to a polynomial-size 
quantum circuilH, [1I1[57], estimating quadratically signed weight enumerators [22j . and estimating average 
fidelity decay of quantum maps [531 
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A Jones Polynomials by Fibonacci Representation 

For any braid b e _B„ we will define Tr(&) by: 

s 



We will use | to denote a strand and \ to denote multiple strands of a braid (in this case n). Qn+i is the 
set of all strings of n + 1 p and * symbols which start with * and contain no two * symbols in a row. The 
symbol 

s 
b 



s 

denotes the s, s matrix element of the Fibonacci representation of braid b. The weight Ws is 

^ _ ( (j) if s ends with p 

^ \ 1 if s ends with *. ^ ' 

(j) is the golden ratio (1 + \/5)/\/2. 

As discussed in [3], the Jones polynomial of the trace closure of a braid b is given by 

Vbtr{A-^) = (-A)3'"(''")D"-iTr(pA(6"')). (11) 

6*'' is the link obtained by taking the trace closure of braid 6. is denotes the writhe of the link 6*''. For 

an oriented link, one assigns a value of +1 to each crossing of the form , and the value —1 to each crossing 
of the form . The writhe of a link is defined to be the sum of these values over all crossings. D is defined 



This includes estimating the trace of the unitary as a special case. 
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hy D = —A? — A ^. PA '■ Bn — > TL„(_D) is a representation from the braid group to the Temperley-Lieb 
algebra with parameter D. SpecificaUy, 



PA{(Ti) ^AE,+A-H 

where Ei . . . En are the generators of TL„(Z?), which satisfy the following relations. 

E,Ej = EjEi for|i-j|>l 
EiEi±iEi — Ei 
Ef = DE, 

The Markov trace on TL„(_D) is a linear map Tr : TL„(_D) C which satisfies 



Tr(l) 
Tr(Xr) 

Tr(X£;„_i) 



= 1 



= Tr(yX) 
= ;^Tr(X') 



(12) 



(13) 
(14) 
(15) 



(16) 
(17) 

(18) 



On the left hand side of equation [TSl the trace is on TL„(Z?), and X is an element of TL„(_D) not containing 
En-i- On the right hand side of equation [TSl the trace is on TL„_i(Z?), and X' is the element of TL„_i(_D) 
which corresponds to X in the obvious way since X does not contain En^i. 

We'll show that the Fibonacci representation satisfies the properties implied by equations [T^ [121 [131 and 
[151 We'll also show that Tr on the Fibonacci representation satisfies the properties corresponding to lTBl [T71 
and 1181 It was shown in 3J that properties I16[ I17i and llSi along with linearity, uniquely determine the map 
Tr. It wiU thus follow that Ti{pP{b)) = Tr^pAib)), which proves that the Jones polynomial is obtained from 
the trace Tr of the Fibonacci representation after multiplying by the appropriate powers of D and {~A) as 
shown in equation 1111 Since these powers are trivial to compute, the problem of approximating the Jones 



polynomial a.t A = e 



reduces to the problem of computing this trace. 



Tr is equal to the ordinary matrix trace on the subspace of strings ending in * plus ^ times the matrix 
trace on the subspace of ^rings ending in p. Thus the fact that the matrix trace satisfies property [T71 
immediately implies that Tr does too. Furthermore, since the dimensions of these subspaces are /n-i and 
/„ respectively, we see from equation that Tr(l) = 1. To address propertv [T51 we'll first show that 



Tr 



V 



\ 



rTr 



(19) 



for some constant 6 which we will calculate. We will then use equation [T^ to relate (5 to D. 
Using the definition of Tr we obtain 



Tr 



b 



V 



) 



1 



/ri0 + /«- 



E 



E 



s p p 



s p * p 

ill 



seQ„-2 I ( \ 
s p * p 



s * p * 

i-L 



E 



spp* s*p* 



s p p p 



FX 

s p p p 



s * p p 



s * p p 



where Qn-2 the set of length n — 2 strings of 
no two * symbols in a row. 



and p symbols which begin with end with p, and have 
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Next we expand according to the braiding rules described in equations [3] and 21 



1 



fn4> + fn- 



E 



E 



/ s * p 



I s p * p 



(j)C 



s p p p 



(f>e 



s p * p 



s p p * \ 



s p p p 



s p p * 



V 



s * p p 



S H: p H= 



S * p p 



We know that matrix elements in which differing string symbols are separated by unbraided strands will be 
zero. To obtain the preceding expression we have omitted such terms. Simplifying yields 



1 



fn + 4>fn- 



I S p 
/ ill 



E 



? P , P 
i 



+ ((/)e + a) 



V 



s p * 



s p p 



s * p 



E (^+'^«) 



s * p 



By the definitions of A, a, 6, and e, given in equation [5l we see that 0e + a = b + (f>a. Thus the above 
expression simplifies to 



/n0 + fn-1 

Now we just need to show that 
and 



E 



s* 



b 

T 

s* 



1 



E (-^6 + a) 



I 



seQ„-i 



1 



b 

T 

sp 



+ fn-1 5 fn-l4> + fn-2 



De + a 



1 



fn4>+fn-l S fn-l(j} + fn-2 

The Fibonacci numbers have the property 

fn4>+ fn-1 _ ^ 
/„_!</)+ /„_2 

for all n. Thus equations and are equivalent to 

1 , 



(20) 
(21) 



and 



tie + a ■ 



(22) 
(23) 



respectively. For A — e '^jt/s ^j^ggg both yield 5 = A — 1. Hence 



Tr 



1 



1 



S fn-l(t> + /ri- 



S* 



E 



^ S* 



E 



I 



b 

T 

sp 
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thus confirming equation 1191 

Now we calculate D from 6. Solving [T^] for Ei yields 



Substituting this into [T51 yields 



(24) 



Tt{X{A-^Pa{<j, 



A-H)) = -Tt{X) 



^ A-'TT{XpA{aij) - A-^TiiX) 
Comparison to our relation Tr(XpA((T,)) = }^^{X) yields 



D 



Tr(X). 



A- 



1 

D' 



^ S 

Solving for D and substituting in A = p-^'^'^l^ yields 

D = 

This is also equal to —A^ — A^^ consistent with the usage elsewhere. 

Thus we have shown that Tr has all the necessary properties. We will next show that the image of the 
representation p_F of the braid group i3„ also forms a representation of the Temperley-Lieb algebra TL„(D). 
Specifically, Ei is represented by 

E,^A-'p^^\a,)-A-H. (25) 



To show that this is a representation of TLn{D) we must show that the matrices described in equation 
[55] satisfy the properties [T31 El and [151 By the theorem of [3] which shows that a Markov trace on any 
representation of the Temperley-Lieb algebra yields the Jones polynomial, it will follow that the trace of the 
Fibonacci representation yields the Jones polynomial. 

Since pp is a representation of the braid group and (TiCTj = (Tjai for — > 1, it immediately follows that 
the matrices described in equation [551 satisfy condition [T31 Next, we'll consider condition [T5l By inspection 
of the Fibonacci representation as given by equation [H we see that by appropriately ordering the basi^ we 
can bring pAicri) into block diagonal form, where each block is one of the following 1 x 1 or 2 x 2 possibilities. 



[b] 



c d 
d e 



Thus, by equation [551 it suffices to show that 



A- 



c d 
d e 



A- 



1 
1 



and 



{A~^a-A 



{A-^b-A 



-2\2 



-2\2 



- D{A-^a- A^"^), 

- D{A-^b- A-^), 



A- 



i.e. each of the blocks square to D times themselves. These properties are confirmed by direct calculation. 



^We will have to choose different orderings for different Ci's. 
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Now all that remains is to check that the correspondence [55] satisfies property [TH Using the rules 
described in equation [4] we have 



b 
a 



e 





d 





a 





d 





c 



*p*p 

*ppp 

*pp* 

pppp 
pp*p 
p*pp 
ppp* 

p*p* 



d 

e 



e d 
d c 







a 



a 
b 



*p*p 

*PPP 

*pp* 

PPPP 
PP*P 
P*PP 
PPP* 
p*p* 



(Here we have considered all four subspaces unlike in equation [6]) Substituting this into equation 1251 yields 
matrices which satisfy condition llSI It follows that equation [25] yields a representation of the Temperley-Lieb 
algebra. This completes the proof that 



for A 



-i37r/5 



B Density of the Fibonacci representation 



In this appendix we will show that pi^\Bn) is a dense subgroup of SU{f„-i) modulo phase, and that 

pl${Bn) and p"p$ (Bn) are dense subgroups of SU{f„) modulo phase. Similar results regarding the path 
model representation of the braid group were proven in [T]. Our proofs will use many of the techniques 
introduced there. 

(4) 

We'll first show that p*J {B4) modulo phase is a dense subgroup of SU{2). We can then use the bridge 
lemma from [1] to extend the result to arbitrary n. 

Proposition 2. p^^} {B4) modulo phase is a dense subgroup of SU{2). 
Proof: Using equation [4] we have: 



b 
a 



*p*p* 
*ppp* 



c d 
d e 



*p*p* 
*PPP* 



We do not care about global phase so we will take 

(") t \ 



(") / \ 
P** {(Ti) 



to project into S'[/(/„_i). Thus we must show the group {A, B) generated by 



A- 



b 
a 



B 



d2 



c d 
d e 



(26) 



is a dense subgroup of SU{2). To do this we will use the well known surjective homomorphism (j) : SU{2) 
S0{3) whose kernel is {±1} (c/. 5 , pg. 276). A general element of SU{2) can be written as 



cos I - I 1 + i sin I - ) [xa^ + yay + za^ 
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where Ox, cry, are the PauU matrices, and x, y, z are real numbers satisfying + + = 1. maps 
this element to the rotation by angle Q about the axis 

X 

V 

z 

Using equations [26l and [5l one finds that and 0(-B) are both rotations by 77r/5. These rotations are 
about different axes which are separated by angle 

6*12 cos-i(2 - %/5) ~ 1.8091137886 . . . 

To show that p'^}(Ba) modulo phase is a dense subgroup of SU(2) it suffices to show that and 0(-B) 

generate a dense subgroup of 50(3). To do this we take advantage of the fact that the finite subgroups of 
50(3) are completely known. 

Theorem 1. (0 pg. 184) Every finite subgroup of SO {3) is one of the following: 
Ck'- the cyclic group of order k 
Dk ■' the dihedral group of order k 
T: the tetrahedral group (order 12) 
O: the octahedral group (order 24) 
I: the icosahedral group (order 60) 

The infinite proper subgroups of 50(3) are all isomorphic to 0(2) or 50(2). Thus, since (t>{A) and 0(5) 
are rotations about different axes, (0(^), <t>{B)) can only be 50(3) or a finite subgroup of 50(3). If we can 
show that {<j){A), 4>{B)) is not contained in any of the finite subgroups of 50(3) then we are done. 

Since 4>{A) and 4>{B) are rotations about different axes we know that {(l){A), 4>{B)) is not Ck or Dk- Next, 
we note that R — (j>{A)^ (j){B)^ is a rotation by 20i2. By direct calculation, 29i2 is not an integer multiple of 
2TT/k for k — 1,2, 3, 4, or 5. Thus R has order greater than 5. As mentioned on pg. 262 of [T7], T, O, and / 
do not have any elements of order greater than 5. Thus, {4>{A), 4>{B)) is not contained in C, O, or /, which 
completes the proof. Alternatively, using more arithmetic and less group theory, we can see that 29i2 is not 
any integer multiple of 27r/fc for any k < 30, thus R cannot be in T, O, or / since its order does not divide 
the order of any of these groups. □ 

(n) 

Next we'll consider for larger n. These will be matrices acting on the strings of length n + 1. These 

(n) 

can be divided into those which end in pp* and those which end in *p*. The space upon which pl^ acts 
can correspondingly be divided into two subspaces which are the span of these two sets of strings. From 
equation |4] we can see that p^Jl\ai) . . . p^Jl\an-3) will leave these subspaces invariant. Thus if we order 
our basis to respect this grouping of strings, pl"''(cri) . . . pl"''(cr„_3) will appear block-diagonal with a block 
corresponding to each of these subspaces. 

The possible prefixes of *p* are all strings of length n — 2 that start with * and end with p. Now consider 
the strings acted upon by pi" These have length n — 1 and must end in *. The possible prefixes of this * 
are all strings of length n — 2 that begin with * and end with p. Thus these are in one to one correspondence 
with the strings acted upon by pi*'' that end in *p*. Furthermore, since the rules 2] depend only on the three 
symbols neighboring a given crossing, the block of p^l\ai) . . . p^J!^ {an-s,) corresponding to the *p* subspace 
is exactly the same as pi" ^''(ci) • . . pi* (ct„_3). By a similar argument, the block of pi"' (cti) . . . pi"' [on-^) 
corresponding to the pp* is exactly the same as pi" ^'(cri) . . . pi* ^'(an-a). 

For any n > 3, pl"'(crn-2) will not leave these subspaces invariant. This is because the crossing cr„_2 
spans the [n — 1)*^ symbol. Thus if the {n — 2)**^ and n^^ symbols are p, then by equation [H pi"' can flip 
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the value of the {n — 1)* symbol. The n symbol is guaranteed to be p, since the (n + 1)* symbol is the 
last one and is therefore * by definition. For any n > 3, the space acted upon by pl"''((T„_i) will include 
some strings in which the (n — 2)*'' symbol is p. 
As an example, for five strands: 



p** (CTi) = 



b 





" 


^p^pp* 





a 





>i=pppp* 








a 


^pp^p* 


a 





" 


*p*pp* 





e 


d 


*pppp* 





d 


c 


*pp*p* 



c d 
d e 
0a 




a 




*p*pp* 
>i=pppp* 
*pp*p* 

*P*PP* 
*PPPP* 
*PP*P* 



We recognize the upper 2x2 blocks of pl*''((Ti), and pfJ{o'2) from equationlHl The lower 1x1 block matches 



(5)/ 



p^^J{(7i) and pfJ{<J2), which are both easily calculated to be [a]. pfJic^s) mixes these two subspaces 



.(3) 



(5), 



We can now use the preceding observations about the recursive structure of {p. 



4,5,6,7...} to 



show inductively that pl"^ (Bn) forms a dense subgroup of SU{fn-i) for all n. To perform the induction step 
we use the bridge lemma and decoupling lemma from 1 . 

Lemma 1 (Bridge Lemma). Let C = A(S B where A and B are vector spaces with dinii? > dim A > 1. Let 
W £ SU{C) be a linear transformation which mixes the subspaces A and B. Then the group generated by 
SU(A), SU{B), and W is dense in SU{C). 

Lemma 2 (Decoupling Lemma). Let G be an infinite discrete group, and let A and B be two vector spaces 
with dim(A) ^ dim(i3). Let pa '■ G SU{A) and pt : G ^ SU{B) be homomorphisms such that Pa{G) is 
dense in SU{A) and Pb{G) is dense in SU{B). Then for any Ua G SU{A) there exist a series of G-elements 
an such that lim„^oo /Oo(Q;n) = Ua and lim„^oo Pb(Q^n) = 1- Similarly, for any Ub € SU{B), there exists a 
series /3„ G G such that lim„^oo Pa{Pn) = 1 and hm„_,oc> PaiM = Ub- 

With these in hand we can prove the main proposition of this appendix. 
Proposition 3. For any n > 3, P^JI\b„) modulo phase is a dense subgroup of SU{fn-i)- 

Proof: As mentioned previously, the proof will be inductive. The base cases are n — 2> and n = 4. As 
mentioned previously, pl*''((Ti) = pl*''(cr2) — [a\. Trivially, these generate a dense subgroup of (indeed, all 
of) SU{1) = {1} modulo phase. By proposition [2l p\^}{ai), and p\^} {(J2) generate a dense subgroup of 
SU{2) modulo phase. Now for induction assume that pi" ^■'(_B„_i) is a dense subgroup of SU{f n-2) and 
pi" ^''(Bn-2) is a dense subgroup of SU{fn-z)- As noted above, these correspond to the upper and lower 



blocks of pl"''((Ti) . . . pi*' ^cr„_2 j. itius, by tne decouplmg lemma, p. 
close to [/© 1 for any U € SU{fn~2) and an element arbitrarily close to 1 © [/ for any U G SU{fn~3)- Since, 
as observed above, pi^\<Jn-i) mixes these two subspaces, the bridge lemma shows that p^^J {Bn) is dense in 
5[/(/„_i). □ 

in] in) 

From this, the density of pip and ppj easily follow. 
Corollary 1. plp^(B„) and Ppl\Bn) are dense subgroups of SU{ fn) modulo phase. 



*■"■* (cr„_2). Thus, by the decoupling lemma, pi^\Bn) contains an element arbitrarily 



Proof: It is not hard to see that 



in)/ \ ln+1) / N 

P*p (<yi) = P** (cti) 



(n) / \ (n+1)/ N 

P*p (o-ri-l) = P** (fn-l) 
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As we saw in the proof of proposition [3l pl"''"^'' (cr„) is not necessary to obtain density in SU{fn), that 
is, (pl"^^''((Ti), . . . , pl"^^''(o-„_i)) is a dense subgroup of SU{fn) modulo phase. Thus, the density of plp^ 
in SU{fn) foUows immediately from the proof of proposition [3l By symmetry, Pp"''(-B„) is isomorphic to 
p*p{Bn), thus this is a dense subgroup of SU{ fn) modulo phase as well. □ 

C Fibonacci and Path Model Representations 

For any braid group -B„, and any root of unity e'^'^/'"', the path model representation is a homomorphism 
from Bn to to a set of linear operators. The vector space that these linear operators act is the space of 
formal linear combinations of n step paths on the rungs of a ladder of height k — 1 that start on the bottom 
rung. As an example, all the paths for n = A, k = 5 are shown in below. 




Thus, the n = A, k = 5 path model representation is on a five dimensional vector space. For k = 5 we 
can make a bijectivc correspondence between the allowed paths of n steps and the set of strings of p and * 
symbols of length n + 1 which start with * and have no to * symbols in a row. To do this, simply label the 
rungs from top to bottom as *, p, p, *, and directly read off the symbol string for each path as shown below. 




* P P P ^ 



s p p *p 



* P *P p 



*PPPP 



In [3], it is explained in detail for any given braid how to calculate the corresponding linear transformation 
on paths. Using the correspondence described above, one finds that the path model representation for k = 5 
is equal to the —1 times Fibonacci representation described in this paper. This sign diff'erence is a minor 
detail which arises only because [3] chooses a different fourth root of t for A than we do. This sign difference 
is automatically compensated for in the factor of (_^)3 writho^ ^^iq,^ both methods yield the correct Jones 
polynomial. 



D Unitaries on Logarithmically Many Strands 

In this appendix we'll prove the following proposition. 

Proposition 1. Given any pair of elements U^p G SU{fk+i) and C/** G SU{fk), and any real parameter e, 
one can in polynomial time find a braid b G Bk with poly(n, log(l/e)) crossings whose Fibonacci representa- 
tion satisfies \\p*p{b) — U*p\\ < e and — C/**|| < e, provided that k = O(logn). By symmetry, the same 
holds when considering pp^ rather than p^,p. 

To do so, we'll use a recursive construction. Suppose that we already know how to achieve proposition [1] 
on n symbols, and we wish to extend this to ti + 1 symbols. Using the construction for n symbols we can 
efficiently obtain a unitary of the form 



Mn-i{A,B) 



A 



B] 



* . 

* . 

* . 

* . 

* . 



PP* 
*PP 
PPP 
P*P 



(27) 
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where A and B are arbitrary unitarics of the appropriate dimension. The elementary crossing cr„ on the last 
two strands has the representatior0 



As a special case of equation [27l we can obtain 



d 
C 



. *p* 
•PP* 
•*PP 
■PPP 
•P*P 



A^diag(a) 



V2 



V2 



V2 



. *p* 
■PP* 
■*PP 
•PPP 
■P*P 



Where < a < 27r. We'll now show the following. 
Lemma 3. For any element 



V21 V22 



one can find some product P of 0(1) Mdiag matrices and Mn matrices such that for some phases (pi and <f>2, 



P 



Vii V12 
V21 V22 



PP* 
*pp 
PPP 
p*p 



Proof: Let i3diag(a) and _B„ be the following 2x2 matrices 



Sdiag(a) 



gja/2 







g-m/2 



and Bn = 



We wish to show that we can approximate an arbitrary element of SU(2) as a product of 0(1) -Bdiag and _B„ 
matrices. To do this, we will use the well known homomorphism <f> : SU{2) 5*0(3) whose kernel is {±1} 
(see appendix [Bjl . To obtain an arbitrary element V of SU{2) modulo phase it sufhces to show that the we 
can use 4>{Bn) and (/)(-Bdiag('3)) to obtain an arbitrary 50(3) rotation. In appendix IB] wc showed that 



a 
b 



and 



e d 
d c 



rotates 



correspond to two rotations of Ttt/S about axes which are separated by an angle of 6*12 ~ 1.8091137886 

By the definition of 0, 0(i?diag(a)) is a rotation by angle a about the same axis that ^ 6 

about. <p{B^) is a TT rotation. Hence, R{a) = (t){Bf^Bdia,g{a) Bf^) is a rotation by angle a about an axis which 
is separated by angleQ 26*12 — tt from the axis that 0(i?diag(a)) rotates about. Q = i?(7r)0(i?diag(a))^(7r) is a 
rotation by angle a about some axis whose angle of separation from the axis that </)(-Bdiag('3)) rotates about 
is 2(26*12 — tt) — 0.9532. Similarly, by geometric visuaHzation, (^(i3diag(a'))Q'/'(-Sdiag(— a')) is a rotation by a 



*Here and throughout this appendix when we write a scalar o in a block of the matrix we really mean al where / is the 
identity operator of appropriate dimension. 

^We subtract tt because the angle between axes of rotation is only defined modulo tt. Our convention is that these angles 
are in [0, tt). 
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about an axis whose angle of separation from the axis that Q rotates about is anywhere from to 2 x 0.9532 
depending on the value of a'. Since 2 x 0.9532 > 7r/2, there exists some choice of a' such that this angle of 
separation is 7r/2. Thus, using the rotations we have constructed we can perform Euler rotations to obtain 
an arbitrary rotation. □ 

As a special case of lemma |31 we can obtain, up to global phase, 



Similarly, we can produce M^J^ . Using Mg. 



swap ' 



. *p* 
■PP* 
■*PP 
•PPP 
■P*P 



and equation [57] we can produce the matrix 



Mr 



C 



* . 



PP* 
*PP 
PPP 
P*P 



for any unitary C. We do it as follows. Since C is a normal operator, it can be unitarily diagonalized. That 
is, there exists some unitary U such that UCU^^ ~ D for some diagonal unitary D. Next, note that in 
equation [27] the dimension of B is more than half that of A. Let d — dim(A) — dini(i?), and let Id be the 
identity operator of dimension d. We can easily construct two diagonal unitaries Di and D2 of dimension 
dim(_B) such that {Di Id){Id © D2) = D. As special cases of equation [27] we can obtain 



1 



and 



and 



Mp 



P 



P 



* . 

* . 

* . 

* . 



>l=pH= 
PP* 
*PP 
PPP 
P*P 

*P* 
PP* 
*PP 
PPP 
P*P 

*P* 

PP* 
*PP 
PPP 
P*P 



where P is a permutation matrix that shifts the lowest dim(i?) basis states from the bottom of the block to 
the top of the block. Thus we obtain 



AI2 = J" swap 



MswanA^_D2^4wip — 



>I=P* 

PP* 
*PP 
PPP 
P*P 
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and 



Di®Id 



* . 



PP* 
*PP 
PPP 
P*P 



Thus 



M1M2 



As a special case of equation [27] we can obtain 



D 



* . 

* . 

* . 



PP* 
*PP 
PPP 
P*P 



Mu 



U 



U 



* . 



•PP* 
•*PP 
■PPP 
•P*P 



Thus we obtain Mc by the construction Mc — MijMiM2M^ . By multiplying together Mc and M„_i we 
can control the three blocks independently. For arbitrary unitaries A, B, C of appropriate dimension we can 
obtain 



Macb = 



A 



C 



B] 



As a special case of equation [28l we can obtain 



M, 



unphasc 



Thus, we obtain a clean swap 



M 



clean 



M, 



unphase-^^swap 



PP* 
*PP 
PPP 
P*P 



. *p* 
•PP* 
•*PP 
•PPP 
•P*P 



(28) 



. *p* 
•PP* 
•*PP 
•PPP 
•P*P 



(29) 



We'll now use A^cican and Macb s-s our building blocks to create the maximally general unitary 



Mgen(V^, W) 



V 






w 



* . 

* . 



*p* 

PP* 
*pp 
PPP 
p*p 



(30) 
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For n+1 symbols, the * . . .*pp subspace has dimension fn^s, and the * . . .p*p and * . . . ppp subspaces 
each have dimension fn-2- Thus, in equation I28[ the block C has dimension /„_2 + fn-3 — fn-i, and the 
block B has dimension /„_2- To construct Mgcn{V,W) we will choose a subset of the /„ basis states acted 
upon by the B and C blocks and permute them into the C block. Then using Macb, we'll perform an 
arbitrary unitary on these basis states. At each such step we can act upon a subspace whose dimension is a 
constant fraction of the dimension of the entire /„ dimensional space on which we wish to apply an arbitrary 
unitary. Furthermore, this constant fraction is more than half. Specifically, /n//n-i ~ 1/0 ~ 0.62 for large 
n. We'll show that an arbitrary unitary can be built up as a product of a constant number of unitaries each 
of which act only on half the basis states. Thus our ability to act on approximately 62% of the basis states 
at each step is more than sufficient. 

Before proving this, we'll show how to permute an arbitrary set of basis states into the C block of Macb- 
Just use Mclean to swap the B block into the * . . .ppp subspace of the C block. Then, as a special case of 
equation [28l choose A and B to be the identity, and C to be a permutation which swaps some states between 
the * . . .*pp and * . . . ppp subspaces of the C block. The states which we swap up from the * . . . ppp subspace 
are the ones from B which we wish to move into C. The ones which we swap down from the * . . . * 
subspace are the ones from C which we wish to move into B. This process allows us to swap a maximum 
of fn-3 states between the B block and the C block. Since /„_3 is more than half the dimension of the B 
block, it follows that any desired permutation of states between the B and C blocks can be achieved using 
two repetitions of this process. 

We'll now show the following. 

Lemma 4. Let m by divisible by 4- Any m X m unitary can be obtained as a product of seven unitaries, 
each of which act only on the space spanned by m/2 of the basis states, and leave the rest of the basis states 
undisturbed. 

It will be obvious from the proof that even if the dimension of the matrix is not divisible by four, and 
the fraction of the basis states on which the individual unitaries act is not exactly 1/2 it will still be possible 
to obtain an arbitrary unitary using a constant number of steps independent of m. Therefore, we will not 
explicitly work out this straightforward generalization. 

Proof: In [23' it is shown that for any unitary U, one can always find a series of unitaries L„, . . . , Li which 
each act on only two basis states such that L„ . . . LiL/ is the identity. Thus L„ . . . Li = U~^. It follows that 
any unitary can be obtained as a product of such two level unitaries. The individual matrices Li, . . . , L„ 
each perform a (unitary) row operation on U. The sequence L„ . . . Li reduces U to the identity by a method 
very similar to Gaussian elimination. We will use a very similar construction to prove the present lemma. 
The essential difference is that we must perform the two level unitaries in groups. That is, we choose some 
set of m/2 basis states, perform a series of two level unitaries on them, then choose another set of m/2 basis 
states, perform a series of two level unitaries on them, and so on. After a finite number of such steps (it 
turns out that seven will suffice) we will reduce U to the identity. 

Our two-level unitaries will all be of the same type. We'll fix our attention on two entries in U taken 
from a particular column: Uik and Ujk. We wish to perform a unitary row operation, i.e. left multiply by a 
two level unitary, to set Ujk = 0. If Uik and Ujk are not both zero, then the two-level unitary which acts on 
the rows i and j according to 



1 



U* U* 
Ujk —Uik 



(31) 



y/\lW+\Uj^\ 

will achieve this. If Uik and Ujk are both zero there is nothing to be done. 

We can now use this two level operation within groups of basis states to eliminate matrix elements of U 
one by one. As in Gaussian elimination, the key is that once you've obtained some zero matrix elements, 
your subsequent row operations must be chosen so that they do not make these nonzero again, undoing your 
previous work. 

As the first step, we'll act on the top m/2 rows in order to reduce the upper-left quadrant of U to upper 
triangular form. We can do this as follows. Consider the first and second entries in the first column. Using 
the operation [31] we can make the second entry zero. Next consider the first and third entries in the first 
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column. By operation [21] we can similarly make the third entry zero. Repeating this procedure, we get all 
of the entries in the top half of the first column to be zero other than the top entry. Next, we perform the 
same procedure on the second column except leaving out the top row. These row operations will not alter 
the first column since the rows being acted upon all have zero in the first column. We can then repeat this 
procedure for each column in the left half of U until the upper-left block is upper triangular. 

We'll now think of U in terms of 16 blocks of size (m/4) x (m/4). In the second step we'll eliminate the 
matrix elements in the third block of the first column. The second step is shown schematically as 




The curly braces indicate the rows to be acted upon, and the unshaded areas represent zero matrix elements. 
This step can be performed very similarly to the first step. The nonzero matrix elements in the bottom part 
of the first column can be eliminated one by one by interacting with the first row. The nonzero matrix 
elements in the bottom part of the second column can then be eliminated one by one by interacting with 
the second row. The first column will be undisturbed by this because the rows being acted upon in this step 
have zero matrix elements in the first column. Similarly acting on the remaining columns yields the desired 
result. 

The next step, as shown below, is nearly identical and can be done the same way. 




I 



The matrix on the right hand side of [32] is unitary. It follows that it must be of the form 



where the upper-leftmost block is a diagonal unitary. Wc can next apply the same sorts of steps to the lower 
3x3 blocks, as illustrated below. 



By unitarity the resulting matrix is actually of the form 



where the lower-right quadrant is an (to/2) x (m/2) unitary matrix, and the upper-left quadrant is an 
(to/2) X (to/2) diagonal unitary matrix. We can now apply the inverse of the upper-left quadrant to the top 
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to/2 rows and then apply the inverse of the lower-right quadrant to the bottom m/2 rows. This results in 
the identity matrix, and we are done. In total we have used seven steps. □ 

Examining the preceding construction, we can see the recursive step uses a constant number of the M„_i 
operators from the next lower level of recursion, plus a constant number of M„ operators. Thus, the number 
of crossings in the braid grows only exponentially in the recursion depth. Since each recursion adds one more 
symbol, we see that to construct MgcniV,W) on logarithmically many symbols requires only polynomially 
many crossings in the corresponding braid. 

The main remaining task is to work out the base case on which the recursion rests. Since the base case is 
for a fixed set of generators on a fixed number of symbols, we can simply use the Solovay-Kitaev theorem |20j. 

Theorem 2 (Solovay-Kitaev). Suppose matrices Ui, . . . ,Ur generate a dense subgroup in SU{d). Then, 
given a desired unitary U G SU{d), and a precision parameter 6 > 0, there is an algorithm to find a product 
V of Ui, . . . ,Ur and their inverses such that \\V ~ U\\ < S. The length of the product and the runtime of the 
algorithm are both polynomial in log(l/(5). 

Because the total complexity of the process is polynomial, it is only necessary to implement the base case 
to polynomially small S in order for the final unitary Mgen(V^, W) to have polynomial precision. This follows 
from simple error propagation. An analogous statement about the precision of gates needed in quantum 
circuits is worked out in ^23] . This completes the proof of proposition [TJ 

E Zeckendorf Representation 

Following [19j . to construct the Fibonacci representation of the braid group, we use strings of p and * symbols 
such that no two * symbols are adjacent. There exists a bijection z between such strings and the integers, 
known as the Zeckendorf representation. Let P„ be the set of all such strings of length n. To construct the 
map z : P„ — !■ {0, 1, . . . , /n+2} we think of * as one and p as zero. Then, for a given string s = s„s„_i . . . si 
we associate the integer 

n 

z{s) = ^Sif,+i, (33) 

where fi is the z**^ Fibonacci number: /i = l,/2 = l,/3 = 2, and so on. In this appendix we'll show the 
following. 

Proposition 4. For any n, the map z : Pn ^ {0, . . . , fn+2} defined by z{s) = X]"=i ^ifi+i bijective. 
Proof: We'll inductively show that the following two statements are true for every n > 2. 

A„ : z maps strings of length n starting with p bijectively to {0, . . . , fn+i — !}• 

B„ : z maps strings of length n starting with * bijectively to {fn+i, ■ ■ ■ , fn+2 — !}• 

Together, An and _B„ imply that z maps P„ bijectively to {0, . . . , /„+2 — !}• As a base case, we can 
look at n = 2. 

pp 
^ 1 
*p «-> 2 

Thus A2 and B2 are true. Now for the induction. Let s„_i £ Pn-i- By equation I33[ 

z[pSn-l) = z(s„_i). 
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Since Sn-i follows a p symbol, it can be any element of Pn-i- By induction, z is bijective on P„-i, thus An 
is true. Similarly, by equation 1331 

z(*S„_i) = fn+l + z{Sn-l)- 

Since s„_i here follows a *, its allowed values are exactly those strings which start with p. By induction, 
An-i tells us that z maps these bijectively to {0, ...,/„ — 1}. Since fn+i + /« — fn+2, this implies i?„ is 
true. Together, An and _B„ for all ti > 2, along with the trivial n = I case, imply proposition [l] □ 
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